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^ : ABSTRACT 
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^ . This paper studies systematically, for the first time, the circumgalactic envi- 

■<sj- . ronment of bright IRAS galaxies as defined by Soifer, et al. (1989). While the 

^ . role of gravitational interaction for luminous and ultraluminous IRAS galaxies 

CN . has been well established by various studies, the situation is by far more obscure 

^ . in the IR luminosity range of the bright IRAS sample, 10^°Lq ^ I^fir ~ 10^^ Lq. 

Q_|. To easily identify nearby companion galaxies, the bright IRAS sample was re- 

Q . stricted to 87 objects with redshift range 0.008 < z < 0.018 and galactic latitude 

+l3 , S ^\ 30° \. A control sample, selected from the Center for Astrophysics redshift 

. survey catalogue, includes 90 objects matching the Bright IRAS sample for dis- 

j> ! tribution of isophotal diameter, redshift, and morphological type. From a search 

^ , of nearby companion galaxies within 250 Kpc on the second-generation Digitized 

Sky Survey (DSS-II), we found that the circumgalactic environment of the Bright 
IRAS galaxies contains more large companions than the galaxies in the optically 
selected control sample, and is similar to that of Seyfert 2 galaxies. We found a 
weak correlation over a wide range of far IR luminosity (10^ Lq ^ LpiR ^ 10^^'^ 
Lq) between projected separation and Lpm which confirms a very close relation- 
ship between star formation rate of a galaxy and the strength of gravitational 
perturbations. We also find that the far IR colors depend on whether a source is 
isolated or interacting. Finally, we discuss the intrinsic difference and evolution 
expectations for the bright IRAS galaxies and the control sample, as well as the 
relationship between starbursting and active galaxies. 



Subject headings: galaxies: statistics - galaxies: interactions - galaxies: star- 
burst - infrared: galaxies 
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1. Introduction 

The IRAS Bright Galaxy Survey (hereafter BIRG survey) by Soifer, et al. (1989) and 
by Sanders et al. (1995) (southern extension) includes all the galaxies brighter than 5.4 Jy 
at 60/xm. The IRAS Bright galaxies are therefore, by definition, the brightest extragalactic 
objects in the sky at dO/im. Prom this survey we learned a wealth of astrophysical informa- 
tion: (1) the far IR (FIR) emission dominates the total luminosity in a significant fraction 
of galaxies; (2) at luminosity log(Ljr/-Lo) ,^11 (the so-called Luminous Infrared Galaxies, 
LIRGs), IR selected galaxies become more numeroTis than optically selected Starburst and 
Seyfert galaxies of comparable bolometric luminosity. At luminosity log{Lir/ Lq) ^ 12, the 
so-called Ultra-Luminous Infrared Galaxies (ULIRGs) exceed the space densities of QSOs 
by a factor of 1.5-2 (Sanders & Mirabel (1996); Sanders, Surace, & Ishida (1999)). 

A considerable number of studies suggest a strong relation between galaxy interactions 
and the highest IR luminosity. ULIRGs are often found to be interacting/merging systems 
(Sanders, Surace, & Ishida (1999)). However, the environment of moderately luminous in- 
frared galaxies (lO^^L© < Lfir < IO^^Lq, hereafter MIRGs), and luminous infrared galaxies 
is not well known yet. In this paper, we study the circumgalactic environment of 87 galaxies 
from the BIRG Survey, with luminosity range 1O^°L0<Lfir< 10^^ Lq. The sample is com- 
posed by MIRGs and a few LIRGs. We also consider whether a correlation may be present 
between FIR properties and the projected separation of the BIRG and its nearest companion 
(§4.6). We then compare the BIRG environment to the one of Seyfert 1 and 2 galaxies (§5.2). 
Finally, we discuss the implications of interaction-induced Lfir enhancement for secular evo- 
lution of galaxies and for the relationship between starbursting and active galaxies (§5). In 
the following discussion, we adopt a Hubble constant of = 75 Km Mpc~^. 



2. Sample Selection 

2.1. Bright IRAS Sample 

The Bright IRAS sample consists of 87 objects, and was compiled from the IRAS Bright 
Galaxy Survey by Soifer, et al. (1989) for the Northern hemisphere, and by Sanders et al. 
(1995) for the Southern one. All objects with galactic latitude > 30° were included in 
the sample. In this way we avoid sampling the galactic plane, where a bias in the detection 
of companions is expected because of both absorption and crowding. We further restrict 
our selection to a volume-limited sample (redshift range 0.008 < z < 0.018). A V/Vmax 
test (Schmidt 1968) gives a value of 0.47 ± 0.05 (rms). Since the BIRG Survey is highly 
complete, this sample is expected to be also complete. 
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The lower z limit (0.008) was chosen to avoid objects with very large angular size, while 
the upper z limit (0.018) was set to include the largest possible number of objects and at 
the same time to avoid very small angular sizes, especially for the companions that could be 
confused with stars (sec § 3.1). It is important to point out that all the objects selected with 
the former restrictions lie in the luminosity range IO^^Lq < Lfir < with MlRGs 

being the wide majority (~ 92% are MIRGs, and 8% are LIRGs). The 60/xm luminosity (in 
ergs A"^) distribution and the LpiR (in solar units) distribution of the BIRG sample are 
shown in the left and right panels of Fig. 1 respectively. 



2.2. Control Samples 

The control sample for this study was randomly extracted from a list of more than 
10,000 objects of the CfA Catalog (Huchra ct al. 1983). This control sample matches: (1) 
the isophotal diameter, (2) the redshift, and (3)the Hubble morphological type distribution 
of the BIRG sample. Only objects with galactic latitude |6j7| > 30° were included. The 
control sample consists of 90 objects. A V/Vmax test (Schmidt 1968) gives a value of 0.48 
d= 0.04 (rms). The CS objects are low infrared emitters as can be seen in Fig. 1. Their 
flux at 60// is usually much smaller than 5.4 Jy, and their luminosity at this wavelength is 
systematically smaller than the luminosity of the BlRGs. Objects without a detection are 
treated as upper limits, using the flux density limits of IRAS. The distribution of upper 
limits is shown by the flUed histograms in Fig. 1. The absolute B magnitude distribution 
was not matched. The B luminosity may be partially correlated with the IR luminosity, since 
both are enhanced through star formation processes. Therefore, any attempt to match the 
B luminosity could bias the control sample towards galaxies with high infrared luminosity, 
which is what we want to avoid. 



3. Analysis 

3.1. Identification of Galcixy Companions 

As in previous environmental studies (Krongold, Dultzin-Hacyan, & Marziani (2001), 
Dultzin-Hacyan et al. (1999a)), the search for galaxy companions was performed automat- 
ically on the DSS-11 with the latest version (1998) of FOCAS (Faint Object Classiflcation 
and Analysis System; Jarvis and Tyson (1981)), and was limited to galaxy companions 
that could be unambiguously distinguished from stars by the FOCAS algorithm. Each set of 
pixels with a flux value larger than the sky threshold is considered an object by FOCAS, and 
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can be classified as galaxy or star only if its diameter is larger than 4 pixels. Since the scale 
of the DSS-II plates is ~ 1.0 arcsec per pixel, the minimum angular size to which FOCAS 
is able to classify objects on the DSS-II is ~ 4 arcsec (which corresponds to ~ 1.4 Kpc of 
projected linear distance). However, we further restrict our search to companion galaxies of 
diameters Dq > 5 Kpc. With our methodology we cannot study smaller objects because 
the distribution of companions is dominated by optical pairs (not physically associated; as 
pointed out in §4.1, optical pairs are the wide majority also in the case of companion di- 
ameter between 5 and 10 Kpc). A third limitation is that FOCAS classifies bright stars as 
galaxies, since they appear as extended objects due to scintillation effects. To avoid gross 
mis-classifications, we checked by eye on the computer screen each object classified by FO- 
CAS as a galaxy. Furthermore, border-line objects of marginally resolved appearance were 
not taken into account to avoid also second order mis-classifications. Effect of plate quality, 
point spread function, sky background, and of automatic identification and measurement of 
companion and background galaxies have been discussed in Krongold, Dultzin-Hacyan, & 
Marziani (2001). They will not be discussed again here; the same effects are still infiuencing 
the analysis of the DSS-II. 

As is customary in many previous works (e.g. Dultzin-Hacyan et al. (1999b); Krongold, 
Dultzin-Hacyan, & Marziani (2001)), the fraction of objects with "physical" companions fphys 
is taken as the fraction with one or more observed companions fohsi reduced by the fraction 
of galaxies with one or more optical companions (derived from Poisson distribution) , namely 
fphys — fobs ~ fopt- The number of background galaxies expected to follow a Poisson statistics 
has been obtained as described by Krongold, Dultzin-Hacyan, & Marziani (2001). 

4. Results 

4.1. Companions Within 3Ds 

We looked for companions in a circular area with radius equal to 3 times the diameter 
of the central object {3Ds)- Our results are summarized in Table 1. 

Companion diameter lOKpc > Dq > 5Kpc Of 87 BIRGs galaxies, ^ 40% has at least 
one companion within 3Ds, vs. 43 % of the 90 objects of the CS. The expected number 
of optical companions from Poisson statistics is 36 %, and 36 % for the BIRG, and CS, 
respectively. If optical companions are subtracted, fphys is ~ 4%, and 6.5% for the BIRG and 
the CS, respectively. These results show that there is not a significant excess of companions 
between Bright IRAS galaxies with respect to non-active galaxies, if all companion galaxies 
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with 5Kpc ^ Dc ^ lOKpc are taken into account. However, this result should be viewed 
with caution since fopt ^ fphys- A statistical approach is not appropriate in this companion 
size range. Any inter-sample difference can be proved as significant only if fphys is estimated 
from redshift measurements for all companion galaxies. 

Companion diameter Dc > lOKpc Of 87 BIRG galaxies, ~ 58.4% has at least one 
companion of diameter D > lOKpc within a search radius 3 Ds, against ~29 % of the 90 
objects of the CS. The expected number of optical companions from Poisson statistics is 
20%, and 18.4% for the BIRG and CS respectively. If fopt is subtracted, fphys is ~ 38.4% 
and 10.9% for the BIRG and the CS, respectively. These results show an excess of large 
companions {Dc > 10 Kpc) in the Bright IRAS galaxies with respect to non-active galaxies. 
A test gives a confidence level for this result of 99.9%. 

4.2. Cumulative Distribution of the Nearest Companion in the BIRG Sample 

and in the Control Sample 

The search radius in all the cases was taken as 250 Kpc of projected linear distance, be- 
yond which we assumed a "non detection." The left hand side of Fig. 2 presents three panels 
with the cumulative distribution of the nearest companion (without correction for optical 
companions) up to a projected linear distance (dp) of 140 Kpc. The upper panel shows the 
cumulative distribution of companions with diameter in the range 5Kpc ^ D ^ lOKpc, with- 
out subtraction of optical com,panions. The middle panel shows the cumulative distribution 
of companions with diameter D > lOKpc, and the lower panel shows the same distribution 
for companions with D > 20 Kpc. The error bars on the control samples frequencies were 
set with a "bootstrap" technique (Efron and Tibshirani, 1993) by randomly re-sampling 
the control sample galaxies into a large number (3000) of pseudo-control samples (i. e., we 
built 3000 pseudo-control samples of 90 randomly-selected galaxies). The uncertainty on 
the companion frequency was set as equal to twice the standard deviation measured from 
the distribution of 3000 companion frequencies computed for each pseudo control sample. 
Comparing the environments of BIRG galaxies and control sample galaxies, it is found that 
there is a statistically significant excess of bright companions {Dc ^ lOKpc) in the infrared 
emitters. For companion diameters 5Kpc < Dc < lOKpc the samples show no significant 
difference. 
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4.3. Distribution of Objects with a "Physical" Companion 

From the Poisson statistics, we calculated fopt at distances 20 Kpc, 40 Kpc, etc. By sub- 
tracting this number from fobs, we built the distribution of the nearest physical companion. 
The right hand side of Fig. 2 presents three panels with this distribution, up to dp ~ 140 
Kpc. The upper panel shows the distribution of "physical" companions with diameter in the 
range 5kpc < Dc < lOKpc. The middle panel shows the distribution of physical companions 
with diameter Dc > lOKpc. The lower panel shows the distribution of physical companions 
with D > 20Kpc. In the latter case, the surface density of objects above this diameter (20 
Kpc) is very low, and the probability of finding optical companions is negligible. Therefore, 
the cumulative fobs is ~ fphys- The error bars on the CS frequencies were again set with the 
"bootstrap" technique. As before, the results show a statistically significant excess of bright 
"physical" companions {Dc > lOKpc) in the BIRG galaxies. For companion diameters 
< lOKpc, there is no significant difference between the two samples. 

4.4. BIRGs vs. Syls and Sy2s 

In order to study the difference between the environment of BIRGs, Seyfert 1, and 
Seyfert 2 galaxies, we used the data obtained for Seyfert environments by Dultzin-Hacyan 
et al. (1999a). The comparison is straightforward since z range of our BIRG sample, search 
radius and diameter hmits are identical to the ones of the Seyfert 2 sample of Dultzin-Hacyan 
et al. (1999a). The cumulative distribution for the projected distance dp of the first observed 
companion for these objects is presented in Fig. 3. The error bars in Fig. 3 were set with the 
bootstrap technique, and are at a 2a confidence level, as in Dultzin-Hacyan et al. (1999a). 
The lower panel of Fig. 3 shows that there is almost no difference in the distribution of first 
companion distances between the BIRG galaxies and the Sy2s. On the contrary, the upper 
panel shows that there is a statistically significant excess (a test gives a confidence level of 
99%) of companions in the BIRG sample with respect to Syl galaxies. A similar difference 
was found between Syls and Sy2s (Dultzin-Hacyan et al. 1999a). 

4.5. Group Membership 

We searched in the environment of our objects to determine whether they belong to 
an association of galaxies. We considered any object with at least two companions with 
diameter Dc > lOKpc within a circle of radius 200 Kpc as member of a group of galaxies. 
25% of the 87 galaxies from the BIRG sample matched the former criteria (43% of the BIRG 



- 7- 



galaxies with at least one companion with Dc ^ lOKpc). Only 4.3% of the control sample 
objects were members of groups as defined here (10% of the CS galaxies in pairs with Dc ^ 
10 Kpc). The results imply that BIRG galaxies are more frequently found as members of 
groups than IR-low emission galaxies. 

Only 14% of the objects in groups belong to compact groups of the Hickson Catalog 
(Hickson, Kindl, & Auman 1989) (this is ~ 3.5% of the 87 objects of the sample). We 
checked whether other BIRG galaxies matched the Hickson criteria, but could not find any. 
BIRGs appear to be preferentially members of groups, although of groups that are looser 
than Hickson's compact groups. 



4.6. Interaction Strength and Infreired Emission 

FIR Luminosity Is FIR emission directly dependent on - or even proportional to - 
interaction strength? Our BIRG sample spans a limited range in Lfir, lO^'' — 10^^ Lq. In 
addition, several objects have companions whose angular separation is less than half the 
maximum width of the IRAS aperture. This implies that a biased correlation could arise 
just because, in the closest pairs, we are measuring the flux of two galaxies. Indeed, if Lpm 
of all small separation ( ^ 1") systems is treated as an upper limit, there is no signiflcant 
correlation between projected linear separation dp and LpiR (and the correlation is significant 
if upper limits are treated as detections!). 

A significant correlation appears only if a wider range of Lfir is considered. We added 
to the BIRG sample data from three samples for which environmental data are available. 
We did not consider systems with a companion whose diameter was 5Kpc ^ Dq ^ 15 Kpc, 
since Fig. 2 shows that most of them may be optical companions. The samples are: 

• Our control sample. 

• The sample of "very luminous" IR galaxies by Wu, Zou, Xia, & Deng (1998), defined 
as galaxies with log Lir ^ 11.15 in solar units. 

• The sample of LIG and ULIG selected by Sanders, Surace, & Ishida (1999), which turn 
out to be composed of early and late stage mergers. 

Fig. 4 shows LpiR versus dp for the galaxies of the above samples. One has to consider 
three major limits to the data: 
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• For several objects, the aperture of the IRAS detectors was larger than the separa- 
tion between IRAS galaxy and nearest companion, making impossible to exclude a 
contribution of the companion to the measured LpiR. 

• The search radius on the DSS-II was hmited to 250 Kpc. There are some objects 
( "isolated" ) for which there is no companion of diameter larger than 5 Kpc within this 
search radius. 

• For several CS galaxies, only upper limits to the fluxes are set. FIR fluxes were not 
available for 11 of 22 galaxies (either isolated or with companion of ^ 15 Kpc). 

All of these limitations introduce a censoring on our data. We considered an upper 
limit to the LpiR of the 11 CS objects that were not detected. For these objects we take 
the flux density limits in the four IRAS bands as upper limits to the source flux density. 
"Isolated objects" were treated as censored in projected separation of the flrst companion, 
and a lower limit to dp was set at our search radius of 250 Kpc ("isolated" sources are the 
ones labelled with horizontal arrows in the lower right side of Fig. 4). Of course, dp values 
are lower limits to the true linear separation, which would be the most meaningful parameter 
to be correlated. However, the effect of chance projection is to spread horizontally the points 
toward the left in the diagram of Fig. 4, but does not go in the sense of creating a false 
correlation. For small dp objects, the infrared flux is measured for the system, and the flux 
for each galaxy is not available. 

We therefore considered LpiR, for small separation objects as an upper limit. We then 
applied the generalized Spearman rank correlation test with inclusion of censored data. The 
correlation appears to be statistically signiflcant, with a correlation coefficient ^ — 0.4 
considering small separation objects as upper limits. The probability of the correlation 
being a chance correlation is ^ 10~^. A best flt using Schmitt's binning regression method 
yields the relationship logLpm/Lo ^ [-0.83 ±0.21] log[dp/{lKpc)] + [11.95 ±0.33] (see Fig. 
4). This result is confirmed by the presence of an analogous correlation between dp and the 
specific luminosity at 60 /xm (plot not shown). 

FIR colors It is important to compare the IR properties of the galaxies with different 
strengths of interaction (and thus, with different projected separation dp). To allow signif- 
icant results to emerge in spite of the bias introduced by projection effects, we considered 
four interaction classes: (1) mergers (5 Sanders objects -|- mergers from BIRG), (2) strongly 
interacting systems (galaxies with companion closer than 30 Kpc (log dp < 1.5, where dp is in 
Kpc), 5 Sanders objects -|- BIRG), (3) weakly interacting systems (galaxies with a compan- 
ion beyond 30 Kpc (log dp > 1.5)), and (4) isolated objects (objects without a companion 
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within our search radius of 250 Kpc). The objects were spht between BIRG and control 
sample of different interaction classes. 

Fig. 5 shows the F(60/xm)/F(100/xm) vs. F(12 /xm)/F(25 /xm) color-color diagram 
for the four interaction classes. Mergers and strongly interacting systems show higher val- 
ues of F(60//m)/ F(100//m) and lower values of F(12 //m)/F(25 //m) while isolated ob- 
jects show lower vahies of F(60//m)/F(100;um) and higher values of F(12 ^m)/F(25 fim). 
Fig. 5 is divided in three regions. In the first one (F(60;um)/F(100/xm) ^ 0.75 and 
F12/F25 ^ 0.65), almost all objects are mergers and strongly interacting. In the second 
region (F(60/im)/F(100/im) ^ 0.75 and F(12 ^m)/F(25 fim) ^ 0.65), there is an agglomer- 
ation of objects of all interaction classes. However the 3 mergers in this region are near the 
border to the first region, and their IR colors are very close to the values of the first region 
mergers. The third region (F(60//m)/F(100/im) < 0.75 and F(12 /im)/F(25 //m) > 0.65) 
shows only objects with a companion beyond 30 Kpc, and isolated galaxies. 

Overall Properties Table 2 reports average and sample standard deviation values of the 
parameters considered in our analysis (Column 1), for different interaction strength classes. 
Columns 2-5 report sample average and sample standard deviation for CS and BIRG isolated 
galaxies {dp ^ 250 Kpc). Columns 6-9 report values for BIRG and CS weekly interacting 
galaxies with dp ^ 30 Kpc. The next columns list the sample average and standard deviation 
for the BIRG sample for the remaining two interaction classes: strongly interacting, and 
mergers (there are no CS galaxies with dp ^ 30 Kpc nor mergers). The last four rows provide 
standard estimates of star-formation-related parameters: (1) the Star Formation Rate (SFR), 
which was computed from LpiR using the standard relationship SFR 4.5 x 10~^''Lp7j:j er^ss-i 
Mq yr~^ (Kennicutt 1998); (2) hydrogen molecular mass (collected from various sources 
in literature and available for 41 objects); (3) the ratio Lfir/Mhs; (4) the depletion time in 
yr simply defined as the hydrogen molecular gas mass assumed over the SFR, — / 
SFR. 

There is a clear continuity on FIR properties and SFR from isolated objects to merg- 
ers (except for the 3 isolated BIRGs, but see below). Lpm increases with the interaction 
strength as indicated by our correlation analysis. Systematic differences in FIR color are 
also appreciable. The depletion time is ^ 10^ yr for a/linteraction classes (including isolated 
objects) in the BIRG sample. In the CS objects, ~ 10^*^ yr, comparable to the Hubble 
time. There is a monotonic trend from isolated galaxies to mergers, in terms of increasing 
SFR and decreasing th2, but it is noteworthy that is not statistically different in the 
various interaction classes. 

Isolated objects from the control sample, and isolated objects from the BIRG sample. 
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have impressively different LpiR. This apparent contradiction needs an explanation. There 
are only 3 isolated BIRGs. NGC 5937, NGC 7083, and NGC 5936 did not show a companion 
larger than 5 Kpc on the DSS-II within 250 Kpc. However, all of these galaxies present 
peculiarities. (1) NGC 5937 has a distorted morphology, and it may have a loop of gas 
which could be a signature of interaction. (2) NGG 7083 is a barred Sc galaxy that hosts 
a LINER. It looks perturbed because of an off-centered loop. (3) NGC 5936 has a highly 
distorted morphology, which may be indicative of recent interaction. These galaxies may 
have been disturbed by the presence of a small companion disrupted or projected over the 
main galaxy. Isolated CS galaxies do not show distortions or peculiarities that could make 
them special objects in terms of morphology or interaction. 

5. Discussion 

The percentage of companion galaxies within 3D5 and the distributions of observed 
and physical companions show an highly significant excess for the BIRGs. The difference 
between BIRGs and CS galaxies is especially striking if large companions with Dc ^ 20 
Kpc are considered (the BIRG galaxies have 3-4 times more companions within fti 140 Kpc; 
strongly interacting systems in the CS may be ^ 1%). Our results also indicate a direct 
relationship between interaction and enhancement of IR emission. We have considered a 
very large range in LpiR, ~ 10^'^ — 10^^'^ Lq, which is unprecedented and probably sufficient 
to overcome the bias introduced by random projection of separation. This may explain why, 
with some notable exceptions (e. g. Sanders & Mirabel (1996) and references therein), 
several previous analyses did not found any convincing correlation between dp and LpiR 
among interacting galaxies. Our result extends to a lower LpiR range and quantifies results 
that were known qualitatively for LIRGs and ULIRGs (Sanders, Surace, & Ishida (1999)). 

5.1. Implications for Star Formation 

An increase in Lpm can be observed across a sequence from isolated galaxies to strongly 
interacting systems. Color variations are consistent with the emergence of a FIR continuum 
component whose luminosity and colors are correlated. This component can be associated 
to thermal re-radiation by dust of hot stars continuum emission. In the most extreme cases 
of isolated CS galaxies, we may have only a cold cirrus component, T~ 20° K. At the 
other end of the FIR color-color diagram, a "hot" component peaking around 100-60 /um 
may have become prominent. The increase in LpiR can be largely ascribed to an increase 
in the SFR, as shown in many previous studies (Kennicutt (1998), and references therein; 
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Sauvage & Thuan (1992)). For the "hottest" sources (F(25/xm)/F(60/xm) > 0.2), however, 
the reprocessed continuum may be due to a non-thermal source (de Grijp et al. 1992). 

The difference in LpiR and Lfir/Mhj (a factor of more than 100 from mergers to isolated 
CS objects, see Table 2) suggests that strong interactions (dp ^ 30 Kpc) are a necessary and 
sufficient condition for an extreme SFR and for a "starburst" (defined as star formation that 
cannot be maintained over the Hubble time), at least for the galaxies of our sample (this 
result may not be generally true if not all mergers of gas rich-galaxies are infrared luminous) 
. A companion that has approached to less than 30 Kpc to a galaxy may need a time 
^ 3 X ^^^d^oKpcV^QQj^^^^i yr to move beyond this distance. The mean depletion time for 
strongly interacting galaxies is ~ 5 x 10^ yr (Table 2). In this case, the interaction time and 
the th are comparable. This means that a galaxy may exhaust its gas before an interaction 
episode is over, on a time much less than the Hubble time. 

On the other hand, the SFRs of weakly interacting galaxies (CS galaxies with dp ^ 30 
Kpc) do not show values that may be considered extraordinary (SFR 0.52MQ/yr). For 
objects whose companion is separated by dp ^ 30 Kpc, the average dp is approximately 112 
Kpc and 67 Kpc in the CS and in the BIRG sample respectively. The SFR is 10 times 
larger in the BIRG than in the CS. This is consistent with tidal forces (ex dp^) driving the 
SFR increase. A weak interaction may produce a moderate enhancement of the SFR of a 
galaxy, but not lead to dramatic effects on its secular evolution. An important implication 
of our results is that at least part of the large dispersion (a factor ~10) for the SFR in 
galaxies of a particular morphological type (sec Kcnnicutt (1998)) may be explained by 
weak interactions (cf Hernandez Toledo, Dultzin-Hacyan, & Sulentic (2001)). 



5.2. Relationship between Star Forming and Seyfert Galaxies 

Our work and many previous ones leave no doubt that gravitational interaction leads 
to an increase of the SFR in gas-rich galaxies. Less clear is the relationship between inter- 
action and the occurrence of non-thermal nuclear activity. In the simplest scheme of Seyfert 
unification, Seyfert 1 and Seyfert 2 are different because of orientation (see e.g., Antonucci 
(1993) for a review): a molecular torus makes obscuration a major factor in the appearance 
of an active nucleus. However, interaction may be a factor leading to the formation of the 
obscuring torus itself, and to the production of extensive circum-nuclear star formation. A 
significant role of interaction introduces an additional degree of freedom (Dultzin-Hacyan et 
al. 1999a) related to environment and, in a broad sense, evolution. 
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5.2.1. The Environment of Seyfert Galaxies 

The main question is then, what is the environment of Seyfert Galaxies? The most 
recent works have found a positive excess of large companions among Seyfert 2s (Sy2s), but 
not among Seyfert Is (Syls) Dultzin-Hacyan et al. (1999a); Laurikainen & Salo (1995); 
de Robertis, Yee, & Hayhoe (1998)). This challenges previous results suggesting an excess 
without differences between Syls and Sy2s (Dahari (1984); RafaneUi, Violate, & Baruffolo 
(1995)). Problems here may arise because of intrinsic inhomogeneity in the discovery tech- 
niques of Sy2 galaxies, as discussed explicitly by Marziani (1991). In addition, it has to be 
taken into account that discovery methods for Sy2 galaxies like the UV excess and the FIR 
color are sensitive to enhanced star formation. For instance, Schmitt et al. (2001) selected 
a FIR-flux limited sample on the basis of FIR color. They found that 31%±10% of Syls 
and 28 %±7% of Sy2s have companions (optical + physical) within three diameters. These 
frequencies are very similar to the frequency found for Sy2s by Dultzin-Hacyan et al. (1999a) 
(companion diameter > lOKpc within 60 Kpc, the case most similar to the one considered 
by Schmitt et al. (2001)). By introducing a bias in favor of star-forming Syls, their selection 
criterion may have simply increased the fraction of interacting Syls (Dultzin-Hacyan et al. 
(1999a) found 21%!) with respect to Seyferts 2. As it can be seen below, almost all AGN 
from the BIRG sample show evidence of significant star formation and belong to interacting 
systems. 

Since an excess is found from uniformly-distributed samples (Dultzin-Hacyan et al. 
1999a) and also for a Seyfert sample selected from the Genter for Astrophysics redshift 
survey (de Robertis, Yee, & Hayhoe (1998)), we consider an excess of bright companion 
among Seyfert 2 galaxies, and no excess among Seyfert Is with respect to a suitably- chosen 
control sample of non-active galaxies as the most accurate representation of the Seyfert 
environment. 



5.3. Star Forming &; Seyfert Galetxies: an Evolutionary Sequence? 

The result of this paper which is relevant at this point is that BIRGs seem to have, to 
a high confidence level, more large and close companions (D > 10 Kpc, d < 60Kpc) than 
Seyfert Is, and seem to be similar in their environment to Seyfert 2s (cf. Dultzin-Hacyan 
et al. (1999a)). This statistical result gives support to a scheme that several workers have 
considered (Heckman et al. (1989), and references therein; Sanders et al. (1988)). The scheme 
is an evolutionary sequence for AGN driven by interaction: 
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Inter action =^ Starburst — > Seyfert 2 — > Seyfert 1 (1) 



where the brace indicates that Syl and Sy2 may be actually the same kind of objects 
seen in different orientations. There are several lines of (admittedly circumstantial) evidence 
that also support this simple evolutionary path: first, the contribution of thermal emission 
to the bolometric luminosity appears to decrease along the sequence (Dultzin-Hacyan & 
Ruano (1996)). Syl nuclei have been revealed in several evolved mergers (for instance, see 
RafaneUi, Marziani, Birkle, & Thiele (1993)). Second, there are several active galaxies in 
the BIRG sample. Of 87 galaxies, 17% host a Seyfert 2 nucleus (15 objects), but only 2.5% 
host a Seyfert 1 nucleus (2 objects). There is no statistical difference between the LpiR of 
active and non-active galaxies, except for a slightly higher value in the Syl objects. The 
value of F(25/im)/F(60/im) for Sy2s is ~ 0.18 and for Syls ~ 0.20, compared to~ 0.13 for 
non-active galaxies. F(25/im)/F(60;um) is larger in Syls and Sy2s due to the contribution 
to the continuum of a non-thermal source (de Grijp et al. (1992)). Estimating in a careful 
way the ratio thermal/nonthermal emission for the BIRGs is not possible from published 
data. However, 73% of the Seyfert 2 (11 of 15) show evidence of significant star formation 
(there is evidence of a circum- nuclear starburst in 45% of the star forming Sy2s). Also the 
2 BIRG Seyfert 1 galaxies show evidence of a circum-nuclear starburst (this makes selecting 
samples of Seyfert 1 and Seyfert 2 from LpiR even more improper for environmental studies 
than selecting them from catalogues!) 

The evolutionary sequence outlined above can be understood in three different ways. 

1. It can be read as a sequence of obscuration properties: (a) fully-obscured Seyfert Is 
(i.e., seen as a Seyfert 2 from all viewing angles), (b) obscuration dependent on viewing 
angle Syls (the "unification" Syl and Sy2 scenario), (c) almost fully unobscured Seyfert 
Is. 

2. It can be a sequence of AGN power; a possibility is that the accretion rate may be 
insufficient to maintain a Broad Line Region (BLR) in some Sy2s. 

3. A low power may also occur in its earlier stages just because the central black hole is 
rather of low mass, maybe because the black hole was not originally present. 

A wealth of X-ray data show that most Seyfert 2 are consistent with an AGN X-ray 
spectrum increasingly less absorbed at energies ^ 5 KeV. This means that an AGN has been 
already switched on in many, if not all, Seyfert 2 galaxies (Moran et al. 2001; Matt 1997). 
This results supports the obscuration sequence (since the power of the AGN will be roughly 
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the same in different types). In this scheme, Seyfert 2s may appear as the low-luminosity 
analogues of ULIRGs, which have been suggested to be precursors of quasars. 

Tran (2001) studied a sample of Seyfert 2 galaxies to determine how many of them 
were obscured Syls i.e., showed a hidden broad hne region (HBLR) in polarized hght. He 
concluded that non-HBLR Sy2s are not more obscured than HBLR Sy2s, but less powerful 
AGN. This result goes against the obscuration scenario, and favors the AGN power scenario 
(points 2 and 3). Gu, Dultzin-Hacyan, & de Diego (2001) studied the properties of 51 Sy2s 
with evidence of high circumnuclear SFR. They found that while Sy2s with a HBLR have 
similar Infrared- Radio properties as Syls, Sy2s without a HBLR have properties similar to 
Starbursts. These results can be straightforwardly understood in the context of an evolu- 
tionary scheme. While objects without a HBLR are "younger" Syls (whether very obscured 
or with very low AGN power), Sy2s with a HBLR are "young" Syls that may keep forming 
stars on their nuclear region but that are less obscured or with higher AGN power. 

Obscuration, low accretion or small black hole mass could be therefore the main physical 
factors behind any evolutionary sequence. However, we think that there is presently not 
enough evidence to decide in favor of one of these factors. 

5.3.1. Environmental Effects as Drivers of any Evolutionary Sequence 

The time needed for Syls to emerge (whether as unobscured or as high power AGN) 
could be longer than the escape time of an unbound companion from the very close environ- 
ment, or comparable to the time-scale needed for an evolved merger (~ 10^ yr). This will 
naturally explain why Starbursts and Sy2s are found more often with closer companions. 

For AGN triggering in a gas rich galaxy, the occurrence of a tidal perturbation may be 
more relevant than its duration (Keel 1996). An hyperbolic encounter may well trigger a 
radial flow in the innermost regions of a gas-rich galaxy. The time needed by the companion 
to move away by 30 Kpc is ~ 1. x lO^dsoKpc^'^^^oxms-i "^^^ timescale for a clump of 
gas to fall from the outer regions of the nucleus (a few hundreds of parsecs) to the inner 
central pc is ^0.1 Gyr (Bekki (2000)), and this can be considered a lower limit to the 
time needed for the onset of the active nucleus. Therefore, an hyperbolic encounter with 
moderate Avr can be such that the companion escape from the close vicinity (^ 60 Kpc) of 
the Seyfert galaxy, leaving a non-interacting Seyfert 1 nucleus. If obscuration is significant, 
or if the AGN power is small (because of low accretion rate, or of an undermassive central 
black hole) , then a longer timescale may be necessary before a Seyfert 1 nucleus is actually 
detected. While BIRGs and Sy2 galaxies have richer environments than Syls at distances 
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^ 60Kpc, the cumulative distribution of the projected separation for the first companion 

(Fig. 3) shows that the environmental difference for Syls, with respect to Sy2s and BIRGs 
decreases dramatically beyond ~ 120 Kpc. This means that, while Sy2s and BIRGs have 
close companions, Syls do have companions, but at higher distances {dp ^ 100 Kpc). Syls 
do not show close companions simply because any activity-triggering interaction took place 
in the past, and, on average, Syl galaxies would not be considered interacting following our 
statistical criterions. 

The limitations of our analysis regarding small companion galaxies (Dc ~ lOKpc) leave 
open other main possibility to account for type-1 activity. It has been proposed that Syls 
may be the result of a "minor merger" which purportedly may not lead to a dramatic star 
formation close to the center of the galaxies and hence to heavy obscuration (de Robertis, 
Yee, & Hayhoe (1998),Taniguchi (1999)). N-body simulations of minor mergers show that 
they produce disturbances in the morphology of the larger galaxy in the first Gyr of the 
onset of the merger, but do not destroy the galactic disk (Walker, Mihos, & Hernquist 
(1996)). Corbin (2000) did not found higher levels of asymmetry in Seyfert galaxies than in 
normal galaxies (in agreement with our work, he found that the most asymmetric galaxies 
were interacting systems with Hll-like spectra). He concluded that, if minor mergers trigger 
AGNs, they appear to do so only in the late stages of the mergers (~ 1 Gyr after the merger 
onset). Minor mergers also boost the star formation of the larger galaxy, but this process 
is not necessarily very dramatic (the induced SFR may be as low as 2MQ/yr) especially 
after the first 0.5 Gyr (Rudnick, Rix, & Kennicutt (2000)). 

The previous mechanisms suggest a revision that complement the unification scheme 
for Seyfert galaxies, and favor the idea of a long timescale to let type-1 AGN emerge. It is 
interesting to stress that times for the onset of this kind of activity are in agreement with 
the time needed to let any unbound companion fly at least few tens of Kpc, or to have a full 
or a "minor" merger (~ 1 Gyr). 

5.3.2. The Unlikely Alternative: No Effect on AGN Triggering by Tidal Forces 

An alternative interpretation for the environmental results for BIRGs, Sy2s, and Syls 
is that interactions may trigger only high SFR but no nuclear, non-thermal activity. Seyfert 
2s may show an interaction-induced enhancement in the SFR as do any other interacting 
galaxies, at least on average. This implies two populations of Sy2 galaxies (Storchi-Bergmann 
et al. 2001): 

• interacting Sy2s with high SFR. The morphology of these galaxies should be distorted 
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due to the interaction, and the interaction could be responsible for the obscuring torus 
of dust, if it exist. The properties of this Sy2 galaxies should be similar to those of 
star forming galaxies. 

• Sy2s isolated and without any circumnuclear starburst. Due to the lack of interactions, 
the morphology of these objects should be very symmetric, without distortions and, 
as Syls, these galaxies should not have any excess of companions when compared to 
normal galaxies. 

If this distinction is correct, strong interactions and non-thermal activity could be fully 
unrelated phenomena (Corbin (2000)). The issue would suffer a 30-year setback. The crucial 
test is then whether the excess of interacting Seyfert 2 with respect to normal galaxies is real 
for a complete sample of Seyfert 2. If not, then there would be no support for a relationship 
between interaction and Seyfert type activity. If yes, then the evolutionary sequence above 
may be appropriate. Defining a complete sample of Seyfert 2 galaxies is tricky, but, as 
already noted, the results based on a limited CfA sample suggest that the difference between 
Seyfert 2 and Seyfert 1 may not be due to sample selection biases. As noted earlier, the 
value of F(25//m)/F(60/xm) for Sy2s is ~ 0.18 and for Syls ~ 0.20. This result argues against 
two Sy2 populations, since Sy2s interacting and with high SFR (like the Sy2s of the BIRG 
sample) should have properties more similar to star forming galaxies, rather than to to Syls 
(and thus, lower F(25/i.m)/F(60/i.m) ratios). 



6. Conclusions 

We studied the environment of Bright IRAS galaxies and compared it to the one of FIR 
low emitters, as well as to the one of Seyfert 1 and 2 galaxies. We found that, on average, 
BIRGs galaxies are more often in interaction, and that their "interaction strength" is higher 
than in a sample of optically-selected galaxies. Our results show a weak anticorrelation be- 
tween the projected separation of the first companion and the FIR luminosity of a galaxy, 
which means an anticorrelation between dp and the SFR. This extends previous results for 
luminous and ultra-luminous FIR galaxies. The FIR properties show a clear and smooth 
continuity as a function of interaction strength, going from very low FIR activity in isolated 
normal galaxies to very high activity in mergers. A consequence is that the FIR luminosity 
function as a function of morphological types is meaningful only for strictly isolated, unper- 
turbed systems. The similar environment found for Seyfert 2 and BIRG galaxies supports 
the possibility of an evolutionary link between Starburst, Seyfert 2 and Seyfert 1. 
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Table 1. Fraction of observed, optical, and physical companions 



Sample Id. 


Sample Size Frequency of Companions (%) 


Significance^ 




Observed Expected Physical 


% 




Companion Diameter > 5 Kpc 


BIRGs 
CS 


87 40.3% 36.3% 4% 
90 42.6% 36.1% 6.5% 


not signif. 


Companion Diameter > 10 Kpc 


BIRGs 
CS 


87 58.4% 20% 38.4% 
90 29% 18.4% 10.9% 


99.9% 



^Statistical significance for the hypothesis that the listed samples are dif- 
ferent from the BIRG sample. 



Table 2. FIR properties for galaxies with different interaction strengths 



Parameter 




ISOLATED 






Separation 


> ZOKpc 




Separation 


< 30i^pc 


MERGERS 


CS 


St. Dev 


BIRGs 


St. Dev 


CS 


St. Dev 


BIRGs 


St. Dev 


BIRGs 


St. Dev 


MERGERS 


St. Dev 




11 




3 




11 




20 




16 




10 




< dp > (Kpc) 


> 250 




> 250 




111.73 


71.50 


66.94 


35.05 


18.20 


7.42 


1.81 


1.17 


<Lfir> (1O"'L0) 


0.21 


0.13 


6.92 


2.53 


0.92 


0.47 


3.83 


2.30 


16.8 


18.8 


71.6 


103.2 


< ii2M >(10:^"ergs-l_ff2-l) 


0.27 


0.15 


1.35 


0.20 


0.34 


0.21 


1.19 


0.71 


5.16 


8.76 


9.43 


1.99 


< i25M >(10»"erss-l_ff2-l) 


0.29 


0.12 


2.98 


1.71 


0.39 


0.30 


2.61 


1.44 


14.6 


22.2 


48.6 


88.4 




1.02 


0.51 


19.4 


8.70 


3.38 


2.80 


20.2 


10.5 


70.3 


67.7 


275.1 


374.2 


< iioOM >(1030ergs-ii/2:-i) 


3.40 


2.27 


42.3 


13.1 


9.85 


6.84 


37.1 


22.0 


103.2 


79.1 


285.3 


344.4 


< F12/F25 > 


0.94 


0.29 


0.53 


0.21 


0.72 


0.16 


0.48 


0.19 


0.41 


0.13 


0.16 


0.05 


< F60/F100 > 


0.33 


0.09 


0.45 


0.09 


0.34 


0.12 


0.57 


0.10 


0.63 


0.17 


0.87 


0.14 


< M(^H2) > (10^ Mq) 


3.77 


1.68 


2.95 


not avail. 


7.12 


2.11 


6.40 


5.89 


7.55 


6.81 


9.60 


9.04 


< Lfir/M^H-z) > {Lq/Mq) 


0.52 


0.42 


26.30 


not avail. 


0.85 


0.53 


13.56 


12.65 


31.47 


35.94 


118.0 


105.77 


< SFR >^ (Mq yr-i) 


0.36 


0.23 


12.10 


4.42 


1.45 


1.15 


6.71 


4.03 


29.56 


32.95 


125.34 


180.09 


< Th^ >b (109 yr) 


19.6 


17.9 


0.22 


not avail. 


8.39 


5.22 


0.91 


0.86 


0.51 


0.59 


0.09 


0.08 



'^SFR « 4.5 X 10-44L^jfl_g^g^3-i M0 yr"! (Kennicutt 1998) 
^T„ = Mh2 / SFR. 
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Fig. 1. — Luminosity at 60 /xm (in ergs A~^) (left panel) and Lpm (in solar units) (right 
panel) for BIRGs and CS galaxies. The sohd line traces the distribution of the BIRG sample, 
and the dashed one that of the control sample, including detections as well as upper limits. 
The filled area identifies the distribution of objects from the CS whose specific luminosity 
at 60 /im and LpiR are known only as upper limits. 
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Fig. 2. — Left: Cumulative distributions of nearest observed companion to BIRG galaxies 
binned over 20 Kpc, with a projected linear distance limit of 140 Kpc. Right: Distributions 
of "physical" companions (corrected for optical companions with Poisson statistics). The 
upper panels show the distributions for galaxies with diameter 5 Kpc < Dc < 10 Kpc, the 
middle panels show "bright" companion galaxies whose diameters are Dc > 10 Kpc, and the 
lower panels show companions with Dc > 20 Kpc. The sohd hue corresponds to the BIRG 
sample, while the dotted-dashed refer to the control sample. The error bars on the CS are 
at a 2cr confidence level. 
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Fig. 3. — The distributions of the nearest companion with diameter Dc > 10 Kpc, binned 
over 20 Kpc, up to a projected hnear distance of 140 Kpc, for Syl, Sy2, and BIRG galaxies. 
Upper panel: BIRG vs. Syl galaxies. Lower panel: BIRG vs. Sy2 galaxies. The sohd hne 
corresponds to the BIRG sample, while the dashed refers to Syl in the upper panel and to 
Sy2 in the lower one. The error bars are set at a 2a confidence level. 
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Fig. 4. — Far Infrared luminosity Lpi^ vs. projected separation dp for Wu, Zou, Xia, & Deng 
(1998), Sanders, Surace, & Ishida (1999), BIRG, and CS galaxies with bright companion. 
107 objects in total. The solid line corresponds to the best fit. 
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Fig. 5. — Color-color plot for objects in different interaction classes, and splitted between 
BIRGs and CS. The figure is divided in three regions. In the first one (F60/F100 > 0.75 
and F12/F25 ^ 0.65), almost all objects are mergers and strongly interacting. In the second 
region (F60/F100 ^ 0.75 and F12/F25 ^ 0.65), there is an agglomeration of objects of all 
interaction classes. The third region (F60/F100 ^ 0.75 and F12/F25 ^ 0.65) shows only 
objects with a companion beyond 30 Kpc, and isolated galaxies. 



